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Abstract 

Last years, there has been a renewed interest in the utilization of statistical field theory meth- 
ods for the description of systems at equilibrium both in the vicinity and away from critical points, 
in particular in the field of liquid state physics. These works deal in general with homogeneous 
systems, although recently the study of liquids in the vicinity of hard walls has also been considered 
in this way. On the other hand, effective Hamiltonian pertaining to the 4 theory family have been 
written and extensively used for the description of inhomogeneous systems either at the simple in- 
terface between equilibrium phases or for the description of wetting. In the present work, we focus 
on a field theoretical description of the liquid vapor interface of simple fluids. We start from the 
representation of the grand partition function obtained from the Hubbard-Stratonovich transform 
leading to an exact formulation of the problem, namely neither introducing an effective Hamiltonian 
nor associating the field to the one-body density of the liquid. Using as a reference system the hard 
sphere fluid and imposing the coexistence condition, the expansion of the Hamiltonian obtained 
yields a usual 4> 4 theory without unknown parameter. An important point is that the so-called cap- 
illary wave theory appears as an approximation of the one-loop theory in the functional expansion 
of the Hamiltonian, without any need to an underlying phenomenology. 
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1 Introduction 



The structure in the interface between two fluid phases at coexistence plays a central role in nu- 
merous specific situations, as for instance wetting transitions and related phenomena [1-7], and 
presents a great interest in itself from a theoretical point of view. Interfaces between coexisting 
subcritical phases at equilibrium provides some of the key features of more complex situations. 
For simple fluids the most important characteristic is the presence of thermally activated capillary 
waves [8] and their coupling with bulk fluctuations. Capillary waves are usualy described through 
the introduction of an effective surface Hamiltonian written as a functional of the interface height 
h{s) (s: coordinate parallel to the interface) and a priori built from a phenomenological description 
of the interface. Then this Hamiltonian, H[h(s)], can be treated in the framework of field theoretical 
methods [2,4], or renormalization group (RG) [9] theory. Recently, the hierarchical reference theory 
(HRT) has been generalized to inhomogeneous case in order to deal with wetting transitions [7]. 

At the liquid vapor interface when H[h(s)] is treated at the gaussian level the capillary wave 
theory (CWT) is recovered, leading to the well known l/q 2 behavior for the height-height structure 
factor Saz Az(q)- Using a 3D </> 4 model [5,6,10,11] the theory has been extended in the sense that 
the stochastic variable of the effective Hamiltonian does not coincide any more with the height h(s). 
An alternative way to deduce an effective Hamiltonian from a DFT functional including the effect 
of the local curvature of the interface has been introduced in refs. [12, 13]. Recently [14], a model 
for the density profile based on an extension of a displaced profile type of approximation, where the 
profile is written as a function of (z — h(s)) was considered. In addition, together with the surface 
fluctuations described by the height h(s)), the bulk phase fluctuations are added. 

The purpose of the present work is to provide a simple description of the liquid / vapor interface 
structure of simple fluids in the framework of statistical field theory. We start from the Hubbard 
Stratonovich (KSSHE) [15] transform to get a representation of the grand-potential from which 
we deduce the surface tension and the surface structure factor at the the so-called one loop ap- 
proximation level. This is done via the determination of the eigenvalues and eigenfunctions of 
the second functional derivative of the mean field KSSHE Hamiltonian. Taking into account the 
whole spectrum of eigenstates is of the outmost importance to obtain the correct result. One of the 
salient points of this work is that one recovers both the CWT and its first extension, namely the 
appearance of the bending rigidity factor k and the coupling between surface and bulk correlations 
in S(q), without invoking any ad hoc phenomenology. 

2 KSSHE transform and mean-field approximation 

We briefly recall the important steps leading to the approriate statistical field theory formulation of 
the problem; the interested reader is refered to ref. [15, 16] and [17] for an alternative formulation. 
We consider a simple fluid whose pair interaction potential includes a hard sphere repulsive part 
and a soft attractive part, denoted v(r). Let w(r) = — /3v(r) as usual and suppose that only purely 
attractive potentials are considered (w is a positive definite operator, i.e. w{k) > 0). In the core 
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region (r < a), w can be chosen at will or regularized in order to fulfill a conveniently chosen 
criterion. We work in the grand canonical ensemble, namely at constant chemical potential v = ftp,, 
and we consider the GC partition function related to the grand potential 5 = exp(— /3fi). Using 
the Hubbard Stratonovich transformation one can get S in the form of a functional integral of 
exp(—H [(/>]) over the field <j)(f) which makes the link with statistical field theory, H[<f>] being the 
effective Hamiltonian 

H= \ < <j>\w- x |0>-ln(H ffs [P + </>]), (1) 

with v = v — w(0)/2 and Shs[ v ] is the hard sphere grand partition function at the local chemical 
potential v{r). A Landau- Ginzburg form is obtained from a functional Taylor expansion of \u(Ehs) 
arround some reference chemical potential of the hard sphere fluid, uq. Performing an expansion 
up to order k 2 of the propagator, A~ 1 (fc) = Kq + -ft^fc 2 , we are left with 

H[<t>] = H + J {^j + §V + V{4>) B^j dr (2a) 

n*) - E u ^ n ( 2b ) 

n>3 



U » = - dn{ !! P n S) (" - = -(PPhs) {n) Wo] (2C) 

where in addition, we have neglected the k 2 dependence of the u n . This formulation allows an exact 
mapping between the densities and their correlations with the mean value and the correlations of 
the field [15]. For instance at coexistence, the densities pi and p g of the two phases correspond to 
the two values < <f) >i and <</>>;. The external field B in f2a]) is given by : 



B = w 1 * Av + phs [vo] with Av = v — v§. 

The value of Vq can be chosen in such a way that the (j) 3 term vanishes and the coexistence condition 
between the liquid and vapor phases is B = 0. Kq is related to the deviation from the critical 
temperature: t — (T c — T)/T c cx (— Ko). The interaction part of the Hamiltonian is denoted by 

In the inhomogeneous system, we assume a special realization of the two phases coexistence: we 
impose explicitly the occurence of a bulk liquid and a bulk gas, at densities pi and p g separated by 
a planar surface located at z = z . The inhomogeneous mean field equation, SH [</>]/<5</> = 0, leads 
to the mean field profile <f> c {z) = 4>b tanh(c(z — z j), when V(4>) is restricted to the order (j) 4 . Here, 
c = (— K0/2K2) 1 / 2 , the inverse of the bulk correlation length, is also the inverse of the intrinsic 
interface width and — ► as t 1 / 2 when t — » 0. The mean field surface tension follows from the 
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identity OmfM = H[<j) c (z)] with the result 

f3a 2 lMF = Wi^KzKl/ul) 1 / 2 , 

which is similar to the expression given by Brilliantov [11] with however another definition for the 
reference system. Note that for t — ► 0, Jmf ~ t 3 / 2 . Finaly, for numerical calculations, we took for 
the u n their Carnahan Starling values and the potential is regularized either in a W.C.A. scheme 
or from ref. [22]. 



3 One loop equations 

In order to go beyond the mean-field level we expand the Hamiltonian in the vicinity of 4> c (z). The 
first correction stems from the second order term and we have (zq = 0), 

H[cp = cp c + X ] ^ H[cb c ] + ^J X (l)tf (2) (l,2)x(2)dld2 (3) 

Introducing the function g(z) = (Kq + V"((f> c (z)))/K2, after a Fourier transform parallel to the sur- 
face, the operator is diagonalized in the set of the eigenfunctions <p n solution of the eigenvalues 
equation (e n = K 2 (q 2 + c 2 uj n )) 

K 2 ( K -§^ + 1 2 + V>n{z) = e u tp n (z) (4) 

the solutions of which are known [5, 18]; the spectrum of eigenfunctions includes two bound states, 
ipo = Co/ cosh 2 (cz) and ip\ = C\ sinh(cz)/ cosh 2 (cz), with luq = 0, u>\ = 3 respectively and a 



subset of unbounded states, or continuum spectrum tfk(z), with k = c^/uk — 4 which behave as 
plane waves in the bulk phases, i.e. far from the interface, and are given by 



<p k (z) = Ci(fc)e* fcz 



k 2 3i 

2 k tanh(cz) + 3 (tanh 2 (cz) - 1) 



c 



(5) 



It is important to note that fo{ z ) oc d<f> c {z)/dz. Then the subset of unbounded eigenstates must 
be orthogonalized; using the eigenvalue equation |(4]) and given the form of the <fk{z) (©, we get a 
non trivial dispersion relation leading to the density of states (x — fc/c) 



m = --- 

TT 7TC 



1 



(1+X 2 ) {A + X 2 )\ TT 



L ■ ^ (6) 



where L is the size of the system in the z direction. It is important to note that f_ f(x)dx = -2 
which proved useful in all calculations. Zittartz [18] have already got this density of states but 
without explicitely mentionning the need of the orthogonalization. In [19] a similar kind of dis- 
persion relation was obtained in the modelling of the charge density profile of electrolytes in the 
framework of another field theory. We emphasize that the orthogonalization of the eigenstates 
is of crucial importance for the calculation of gaussian functional integrals. In order to calculate 
ln(H) and the correlation functions, we write the fields %(r) in the basis which diagonalize the 
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operator flw, namely a Fourier transform in the direction s*and a projection on the fx(z) where 
{A} = {n = 0,1; k} denotes the whole spectrum of eigenstates. From usual gaussian functional 
integrals [20] we get 



In(H) 



HI 



d 3 k 
(2tt) 3 



In 



w-\k) 



(7) 



^(E *(*(*))+/ 



dkf(k/c) \n(e k (q)) 



where e& corresponds to the second functional derivative of the Hamiltonian for the 4 homogeneous 
model taken at 4> = 4>c(i L). Therefore, the volume term of ln(S) in is nothing but f3pV at 
the one loop level of approximation [15]. We finally write the surface term of ln(E) in term of the 
surface tension 7 = [ln(S 6 ) - ln(3)]/S = 7a/f + 7 (1) - We obtain (7 (1 * :) = /?7 (1 V 2 ) 



' 8tt 



+ (.T-4c 2 )ln 1 - 



(a; - c 2 )ln ( 1 
4c 



2c 



\/i+2c 



— < oc\fx 



X + c 



(8) 



where q m and qm are the lower and upper bound respectively of the integral over q. A similar 
result was obtained in [18] for a spin model. This result differs from that obtained in [4] where only 
the n = eigenstate is kept. 

We now consider the calculation of the field two-body correlation functions, 
G(l,2) =< x(l)%(2) >h- More precisely, we focus on the Fourier transform parallel to the sur- 
face, G(z\,Z2,q)- For this we have to calculate the sum 



E 



tp X {zi)(p* x {z2) 

£ a(<7) 



E 

71=0,1 



<Pn(zi)<p n (z 2 ) 

£n(q) 



P {k f k{Zl) f\ {Z2) dk 

efe(<?) 



A contribution of the integral over k cancels exactely the direct contribution of the two bound states, 
which shows, once again, that the approximation consisting in keeping only the tpo(z) eigenmode 
is not sufficient. The result is (z>,< = sup,inf(zi, z 2 )) 



G(zi,z 2 ,q) 



9 exp(-c|z 12 |(( g /c) 2 +4) 1 / 2 ) 
2c^ 2 ((q/c) 2 + 4)V2( ff / C )2(((g/c)2 + 3) 

1 + q 2 /3c 2 + ((q/c) 2 +4) 1/2 tanh(cz>) +tanh 2 (c2>) 
1 + q 2 /3c 2 - ((q/c) 2 +4) 1/2 tanh(cz<) + tanh 2 (cz<) 



(9) 



4 Results and discussion 

Capillary behavior and surface structure factor 

We started from the expansion of the effective Hamiltonian on the basis of the eigenstates ip\ (z). 
The first eigenstate, (po(z) is proportional to the derivative of the mean field result 4> c {z), the 
proportionality constant being determined by the normalization. If we keep only this first eigenstate, 
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the expansion of x reads x(s, z) = £(s)tp (z). Hence the field takes the form (f>{z) = <j) c {z — z int (s)) 
and only the linear term in the expansion of 4> with respect to Zi„t is kept. This corresponds to the 
so-called rigidly displaced profile approximation, where Zi nt — — a£(s) represents the fluctuating 
location of the interface. The corresponding contribution to H is a functional of £(s) which defines 
an effective surface Hamiltonian given by 

Hi Q) W)} = \J ds[K 2 (dMs)) 2 +u a m 2 ] (10) 
using uj = 0, and the mean field equation yielding K 2 /a 2 = fi^juF we rewrite iflOj) in the form 

#i 0) = \PlMF J ds{d S {z lnt {s)f = HcWT (11) 

which coincides with the usual effective surface Hamiltonian of the CWT theory for the free surface, 
i.e. in the absence of external field. Therefore we obtain the CWT as the lowest approximation 
beyond the mean field approximation without invoking phenomenological arguments. 

The structure is also characterized by the height-height correlation function, < Zi n t(s-\)zi n t (s^) > 
or its Fourier transform parallel to the surface which defines the surface structure factor, SazAz(q), 
where Zi nt (s) is the location of the interface relative to its mean value. We consider J x(s, z)dz/A(pb, 
where A<f)b = <f> c {L) — <t>c{—L), as a measure of the instantaneous location of the surface at s, which 
amounts to define the location of the surface from a constraint on the integral of x, as is done in [14] 
for the density profile. We are then led to identify Saz Az(q) = (A4>c{b))~ 2 J G(z±, z 2 , q)dzidz 2 , 
which corresponds to the Si C used in [14]. It is important to notice that the coupling with the bulk 
fluctuations are included in the present formulation through the eigenstates of the continuum. We 
also introduce the effective surface width, or surface corrugation, a e f* =< Zi nt (si)zi n t(si) >■ 

The behavior of S&z Az(q) is analysed from the function g(q) = j dz\ J dz 2 G(zi, z 2 ,q). From 
© we get for the leading term of G(z\, z 2 , q) when q — > 

G(zx,Z2,q^> 0) ~ ■^-^(po(z 1 )ip (z 2 ) (12) 
K 2 q z 

which corresponds to the CWT behavior. Then it is easy to show that l/(A(j)b) 2 g(q — * 0) — ► 
1/(7mf« ) from the relation already used between the normalization of ifio(z) and Jmf ((/ <fio(z)dz) 2 
= (A4>i > ) 2 K 2 /jmf)- This is exactly the CWT behavior, leading to the well known logarithmic di- 
vergence of the surface corrugation for which we get a e ^ = (4ttjmf)~ 1 ^(Qm / Qm) ■ 

It can be shown from |(9]) that the contribution to g(q) diverging with the system size, due to 
the bulk correlations, is exactly (2L) times the integral of Gb(zi2, <?),the correlation function of the 
bulk phase, over z\ 2 If we keep only these two terms we get 

2L f 1 

Saz Az(q) ^ -t-t^ / G b (q,zi 2 )dz 12 H (13) 

J "fMFq 

Given that the bulk term vanishes when q — ► 0, we see that 11131) presents a cross-over like behavior 
in terms of wave vector q, where a threshold value q s naturally appears, separating the capillary 
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wave behavior at small values of q from the bulk like behavior at q > q s , with q s given by 

q s ~ [(2i 7Mf /A$ 2 ) J G b (q,z 12 )dz 12 ]- 1/2 . 

Such a behavior is in agreement with that of Refs. [14,21]; however, in the present formulation 
the bulk fluctuations come out of the calculation through the continuum subset of eigenstates and 
have not to be added to the field profile. Furthermore, from l|13p we can drop exactly the bulk 
contribution, and doing this we define a purely interfacial contribution to Saz Az(q) (see fig.JI])). 

S int {q) = S Az Az (q) - (2L/A0 2 ) J Q b (q, z 12 )dz 12 (14) 

Then the departure of Sm (q) from its 1/q 2 behavior allows us to isolate the deviations from the 
capillary wave like behavior of S AzAz (q). To this aim we define (3{q) = (q 2 Smtiq)) -1 —jmf- 
/3(q)/q 2 is found nearly constant and thus lead us to define a bending rigidity of the interface, 

K = lim((l/ g 2 )/3( g )) ~ ((l/q 2 )P(q)) (15) 
<7-»0 

which takes a positive value, as it should be for the stability of the interface [13]. In other words, 
the small q behavior of Si nt {q) is 

Si n t{q) ^ i/Ita/f q 2 + nq 4 ] 

Here we find K = k*(7mfc 2 )/(cc) 2 , with the reduced value n* = 0.288. 

10 5 | 1 1 1 1 1 1 1 1 1 r-| 




1Q -1 I , , , ■ ■ J 

0.01 0.1 1 

Figure 1: Log-log plot of g*(x) versus x = q/c for Lc = 100, 80, 60, 40 
and 20 from top; interfacial contribution A^S* nt (x), bottom; 
CWT limit, straight line. 
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Surface tension 

The contribution 7^ 1 - ) , given by ([8]), depends on the two bounds q m and qM of the wave vector q 
which are related to the relevant parameters of the interface: on the one hand qujc = 2-njiac) 
where l/(cc) is the intrinsic width of the interface in unit cr, qu/c £ [l,oo[. On the other hand, 
<7m/c = (2ir)(L x /c) where L x is the system lateral size; hence q m /c = {qM / c){<j / L x ). Thus qM/c 
appears as a natural parameter, with qujc <= [27r,oo] and \/{qM/c) oc c oc \fi. We can now 
re-write 7^*' is a more convenient form (a/L x £ ]0, 1]) 



7 (i*) 



2( 9 m/c) 2 



7(4 +Wc) 2 )- 7(4 +C—^) 2 ) 

c ±j x 



(16) 



(L x ja) is either the actual lateral system size or the scale at which 7 is measured, for instance 
in a numerical simulation (see ref. [23]), but in any case does not depend on t. We can note that 
whatever the values of a/L x or (qM/c), a 2 "/^ remains finite and more precisely 

which — > when t — > as t 1 / 2 . The results for tr 2 7^ 1 *' are displayed in fig. J2]). We interpret 
o 2 ^ 1 ** 1 {o / L x ) as the q— dependent contribution to 7 with q = (2na/L x ). This means that the 
flat interface corresponding to the mean field approximation is obtained when no fluctuation at all 
are taken into account, namely for q ~ qM- This differs from what is done in Ref. [13] (see also 
ref. [24]) where the contribution to 7 due to the surface fluctuations vanishes at q — > 0. The small 
q behavior of 7^*' is easily obtained from lfT6|) and yields 

7o^o - 7 (1 * } (°) + IT Kw) 2 [a + m ( ccr )] - 2 (9 CT ) 2 Mw)} with a = 3.48491 (18) 

from which we can deduce a crossover value of q given by q^a = (ccr)e a / 2 , separating the q 2 behavior 
from the (qa) 2 ln(gcr) behavior obtained for q < q n and q > qo respectively. It is important to note 
first that we always get an increasing 7^(9) an d secondly that since qo is proportional to c, we get 
a plateau (corresponding to the q 2 dependence) at small values of a/L x only when qM/c takes small 
values, i.e. for the lowest temperatures (see figure ([2]). This behavior is in qualitative agreement 
with the simulation results of [23]. 
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Figure 2: ^{a / 'L x ) /7 (1) (0) in terms of a/L x for q M /{2-Kc) = 1; 2; 5 
and 50 from bottom to top. 

Notice that the term proportional to q 2 in the variation of 7W with q can be interpreted as 
resulting from the energy necessary to bend the interface, and should be related to K obtained from 
the behavior of Si n t (q) (see eq. (fl5|) . This is not a priori the case since we do not expect a fully 
coherence in the framework of a loop expansion between the energetic and structure quantities. 

Concluding remarks 

To conclude, we have shown in this work that the simple one-loop expansion of the grand potential 
in the the simplest inhomogeneous situation provides a qualitative but nevertheless coherent picture 
of the structure of the liquid vapor interface of simple liquids. We emphasize that the spectrum 
of eigenstates resulting from the diagonalization of the second functional derivative of the effective 
Hamiltonian must be treated as a whole and one cannot take into account only the first bounded 
state as sometimes done in the literature. Moreover, starting from a 3D Hamiltonian, and without 
invoking a phenomenological description of the interface, we recover the usual CWT surface Hamil- 
tonian as the first approximate step. We also show that we can extract a rigidity bending factor 
which takes a positive value, in agreement with the requirement for the stability of the interface 
with respect to fluctuations. 
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